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Abstract 

Two-component signal transduction systems, where the phosphorylation state of a regulator protein is modulated by a 
sensor kinase, are common In bacteria and other microbes. In many of these systems, the sensor kinase is bifunctional 
catalyzing both, the phosphorylation and the dephosphorylation of the regulator protein in response to input signals. 
Previous studies have shown that systems with a bifunctional enzyme can adjust the phosphorylation level of the regulator 
protein independently of the total protein concentrations - a property known as concentration robustness. Here, I argue 
that two-component systems with a bifunctional enzyme may also exhibit ultrasensitivity if the input signal reciprocally 
affects multiple activities of the sensor kinase. To this end, I consider the case where an allosteric effector inhibits 
autophosphorylation and, concomitantly, activates the enzyme's phosphatase activity, as observed experimentally in the 
PhoQ/PhoP and NRII/NRI systems. A theoretical analysis reveals two operating regimes under steady state conditions 
depending on the effector affinity: If the affinity is low the system produces a graded response with respect to input signals 
and exhibits stimulus-dependent concentration robustness - consistent with previous experiments. In contrast, a high- 
affinity effector may generate ultrasensitivity by a similar mechanism as phosphorylation-dephosphorylation cycles with 
distinct converter enzymes. The occurrence of ultrasensitivity requires saturation of the sensor kinase's phosphatase activity, 
but is restricted to low effector concentrations, which suggests that this mode of operation might be employed for the 
detection and amplification of low abundant input signals. Interestingly, the same mechanism also applies to covalent 
modification cycles with a bifunctional converter enzyme, which suggests that reciprocal regulation, as a mechanism to 
generate ultrasensitivity, is not restricted to two-component systems, but may apply more generally to bifunctional enzyme 
systems. 
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Introduction 

Two-component systems (TCSs) are modular signal transduc- 
tion systems which are utilized by bacteria and other microbes to 
respond to intracellular or environmental stimuli [1,2]. 'Classical' 
TCSs consist of a sensor histidine kinase (HK) and a cognate 
response regulator (RR), which often acts as a transcription factor 
to activate or repress a particular set of response genes. Upon 
stimulation, the HK autophosphorylates at a conserved histidine 
residue and transfers the phosphoryl group to an aspartate residue 
in the receiver domain of the RR. Often, the unphosphorylated 
form of the HK also exhibits phosphatase activity towards the 
phosphorylated form of the RR (RR-P) endowing many HKs with 
a bifunctional desigTi (Fig. 1). In addition, some RRs exhibit 
intrinsic phosphatase activity which leads to autodephosphoryla- 
tion of RR-P with a half-life ranging between seconds to hours [1] . 

Even though the overall signal flow from the sensor kinase to the 
response regulator is well-conserved between difFerent systems 
there exist substantial variations in the particular mechanism 
through which the phosphoryl group is transferred to the regulator 
protein [3] . To better understand their regulatory properties it has 



become a useful strategy to compare difFerent TCS architectures 
based on their potential input-output behavior. Following that 
strategy, it has been argued that phosphorelay systems, where the 
phosphotransfer to the RR does not occur in a single step but via 
additional intra- or intermolecular reactions [4], may generate 
ultrasensitivity and robustness against noise [5]. Systems with a 
split histidine kinase comprise another class of TCSs where a 
functional HK is generated through binary association between 
two distinct proteins each of which alone would not be able to 
phosphorylate the cognate RR(s) [6] . A theoretical study suggested 
that such systems can potentially exhibit ultrasensitivity and 
bistabUity if the phosphatase activity is predominantly located on 
the free form of one of the proteins making up the split kinase [7] . 
Yet another study compared TCSs with a mono- and a 
bifunctional HK arriving at the conclusion that ultrasensitivity 
and bistabUity can also occur in classical TCSs if the unpho- 
sphorylated forms of the HK and the RR form a dead-end 
complex and if the dephosphorylation of the RR mainly occurs via 
an HK-independent phosphatase [8] . 

In contrast, systems with a bifunctional design are expected to 
generate graded responses to input stimuli [8-10] and to mediate 
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Author Summary 

Bacteria often use two-component systems to sense and 
respond to environmental changes, which involves auto- 
phosphorylation of a sensor kinase and phosphotransfer to 
a cognate response regulator. However, despite conserva- 
tion of this 'classical' scheme there exist substantial 
variations in the mechanism of phosphotransfer among 
systems. Also, many sensor kinases exhibit phosphatase 
activity raising the question whether such a bifunctional 
architecture enables special regulatory properties in the 
response behavior to input signals. According to previous 
studies, classical two-component systems are unlikely to 
produce sigmoidal response curves (ultrasensitivity) if the 
sensor protein is bifunctional. Here, I argue that this is not 
necessarily true if the input stimulus (allosteric effector) 
reciprocally affects multiple activities of the sensor kinase, 
as it seems to be common for bifunctional enzymes. To 
this end, I propose and analyze an extension of the 
experimentally well-supported Batchelor-Goulian model 
which shows that ultrasensitivity requires a high-affinity 
effector and saturation of the phosphatase activity. The 
underlying mechanism involves sequestration of the 
effector by the sensor kinase which restricts the occur- 
rence of ultrasensitivity to sufficiently low effector 
concentrations. Hence, this operating regime might be 
useful to sense effector limitations or to amplify weak 
Input signals. 

concentration robustness [1 1-13]. The latter property means that 
the system response (concentration of phosphorylated RR) is 
invariant with respect to variations of the total RR and HK 
concentrations, at least in a certain range of concentrations. 
Moreover, based on theoretical studies of covalent modification 
cycles with a bifunctional converter enzyme it has been argued 
that ultrasensitivity is unlikely to occur in such systems if the 
bifunctional enzyme employs only a single catalytic site for its 



opposing activities [14,15]. Based on this argument it, thus, 
appears unlikely that classical two-component systems with a 
bifunctional sensor kinase would exhibit ultrasensitivity given that 
the phosphotransferase and phosphatase activities of the sensor 
kinase are believed to occur on a single catalytic site in the 
dimerization domain of the protein [16,17]. Interestingly, this 
conclusion does not apply to bifunctional enzymes with two 
distinct catalytic sites where ultrasensitivity may arise from the 
formation of a ternary complex between the enzyme and its two 
substrates [18] as observed experimentally in the uridylylation 
cycle of the PII protein [19]. 

In the present study, I wish to argue that ultrasensitivity may stiU 
occur in two-component systems with a bifunctional enzyme 
kinase if the input signal reciprocally affects multiple activities of 
the sensor kinase. Reciprocal regulatory patterns have been 
observed in the PhoQ/PhoP system which mediates adaption in 
response to Mg^+ limitation as well as in the NRII/NRI system 
which mediates adaptation to nitrogen limitation by sensing the 
concentration of deuridylylated PII protein in the cytosol. In both 
cases, binding of an allosteric effector (Mg^+ or PII) inhibits the 
autokinase activity and, concomitantly, activates the phosphatase 
activity of the respective sensor protein (Fig. 2A) [20,21]. Indeed, 
based on structural analysis of HK domains it has been argued 
that reciprocal regulation could be quite common for bifunctional 
enzymes [17]. 

In a first step, the impact of reciprocal regulation is analyzed in 
covalent modification cycles with a bifunctional converter enzyme, 
which will serve as a 'toy' model that allows for an intuitive 
understanding of the potential mechanism for the generation of 
ultrasensitivity. In a second step, it wiU be shown that the same 
mechanism may also generate ultrasensitivity in classical TCSs 
with a bifunctional sensor kinase. To this end, an extension of the 
experimentally well-supported Batchelor-Goulian model (see 
below) is proposed which assumes that autokinase and phospha- 
tase activities of the HK are reciprocally regulated by an allosteric 
effector (Fig. 2C). Analysis of this model shows that if the affinity of 
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Figure 1. Signal flow in classical two-component systems. Typically, the sensor histidine kinase (HK) is a (dimeric) transmembrane protein 
which senses extracellular signals directly or through their concentration in the periplasm [3]. In some case, signal-sensing may also occur in the 
cytosol or in the plasma membrane [43]. The HK exhibits up to three distinct activities: (1) autokinase activity leading to the autophosphorylation of 
the HK, (2) phosphotransfer to the response regulator (RR) and (3) phosphatase activity towards the phosphorylated form of the RR {RR' ). In general, 
the input signal may affect all three HK activities although autokinase and phosphatase activities appear to be the most common targets of 
regulation [20,21,44,45]. The phosphorylated form of the response regulator often acts as a transcription factor which activates or represses a 
particular set of response genes including those of the RR and the HK themselves (autoregulation). 
doi:1 0.1 371/journal.pcbi.1 00361 4.g001 
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Figure 2. Reciprocal regulation in two-component systems. (A) Schematic representation of reciprocal regulation in the PhoQ/PhoP [20] and 
NRII/NRI systems [21], In both cases, an allosteric effector (Mg-+ or Pll) inhibits autophosphorylation of the sensor kinase and increases the enzyme's 
phosphatase activity. (B) Batchelor-Goulian model [11] based on the three activities of the sensor kinase (cf. Fig. 1): (1) Autophosphorylation of the 
sensor kinase (HK), (2) phosphotransfer to the response regulator (RR) and (3) dephosphorylation of the RR. Cofactors such as ATP are assumed to be 
constant. (C) Extension of the Batchelor-Goulian model to include reciprocal regulation of the HK's activities as schematized in (A). Binding of the 
allosteric effector L (4) inhibits autophosphorylation (1 ) and activates the phosphatase activity (3) of the sensor kinase. For simplicity, the free form of 
the enzyme (HK) is assumed to have no phosphatase activity vj/hereas the effector-bound form (HKl) is assumed to have no autokinase activity. 
doi:1 0.1 371 /journal.pcbi.1 00361 4.g002 



the efTector is low (as in the case of ]V[g^+ for PhoQ) the system 
exhibits a graded response to changes in the eflfector concentration 
and stimulus-dependent concentration robustness - in agreement 
with experiments in the PhoQ/PhoP system [22]. In contrast, a 
high-affmity effector may lead to ultrasensitivity at low effector 
concentrations, but requires saturation of the sensor kinase's 
phosphatase activity. Comparison of the model predictions with m 
vitro experiments suggests that in the NRII/NRI system the 
occurrence of ultrasensitivity is (partly) suppressed by the intrinsic 
autophosphatase activity of NRI. 

Concentration robustness in the Batchelor-Goulian 
model 

To rationalize the occurrence of concentration robustness in the 
EnvZ/OmpR system of £. cofi, Batchelor and Goulian proposed a 
simple mathematical model based on the three activities of the 
bifunctional EnvZ (denoted by HK in Fig. 2B). Guided by the 
observation that the total OmpR concentration is much larger 
than that of EnvZ [23] ([OmpR].j./[EnvZ].j. w35) they have 
argued that, in the limit [EnvZ].j./[OmpR]y^0, the steady state 
concentration of OmpR-P (denoted by [i^i^^] in Fig. 2B) is 
determined by a quadratic equation [1 1], which can be written in 
the form [SI Text ST) 



RR''f-[RR''] (Rt- 



CpRr '- 



-0. 



:i) 



[OmpR-P] exlribits (concentration) robustness with respect to 
changes in EnvZx. 

Interestingly, Eq. (1) also predicts concentration robustness of 
[i^i^^] with respect to the total concentration of the response 
regulator {Rt) under certain conditions. To see this more 
explicidy, it is worth mentioning that a structurally similar 
equation has been analyzed previously in the context of 
concentration robustness for covalent modification cycles with a 
bifunctional converter enzyme [24]. This analysis has shown that 
the shape of the stimulus-response curve, described by Eq. (1), 
depends on the relative magnitude between the two parameters Cp 
and C, [18]. To this end, it is useful to consider two limiting cases 
corresponding to Ct«Cp and Ci»Cp. In the first case, the 
physiologically reasonable solution of Eq. (1) can be approximated 
by {SI Text SI) 



( 
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cji- 



c, 



Cp — J^T 



Rt < Cp 



c, 



(2) 



Rt-Cp) Rt>Cp 



whereas, in the second case, one obtains the approximate solution 



c,»c, 



p 



[RR'' 



CpRr 



C„ + C, + Rt' 



(3) 



Here, Rt denotes the total OmpR concentration, and the 
parameters Cp = (k^ / kp)Kp and Q = (k^ / kt)K, are proportional 
to the Michaelis-Menten constants associated with the phospha- 



total EnvZ concentration (EnvZx). Hence, the Batchelor-Goulian 
model predicts that, in the limit [EnvZ] j/ [OmpR] y->0, the 
concentration of OmpR-P is approximately independent of 
variations in the total concentration of the sensor kinase, i.e. 



In any case, from the expressions in Eqs. (2) and (3) it is readily 
apparent that [RR''] becomes independent of the total RR 
concentration if the latter is sufficiently large, i.e. [i^i?^] 



tase {Kp = {kp+k2 )/k2) and pho.sphotransferase (Ki = {k,+ RT»Cp (Eq. 
k^)/k^) reactions. Note that Eq. (1) does not depend on the parameter Cp 



2) or Rt»C, (Eq. 3). 
determines both, the 



'Cp if 
Hence, if C,«Cp, the 
threshold concentration 



beyond which [i^i^ ] becomes approximately constant as well as 
the value of that constant. In contrast, if C/»C^, the predicted 
threshold concentration (C/) is much larger than the asymptotic 
phosphorylation level of the response regulator (Cp). Also, the 
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approach to the asymptotic level is different for the two regimes: If 
C(« Cp, [RR^\ increases approximately linearly with Rj up to the 
threshold (Eq. 2) whereas, in the opposite case, it increases 
hyperbohcaUy (Eq. 3). Due to the linear relationship between 
and Rj in Eq. (2) the regime Q « Cp has been called 
'signal-transducing' in Ref [25]. 

Together, Eqs. (2) and (3) suggest that there exist two different 
regimes for the occurrence of concentration robustness and, as will 
be shown below, there is experimental evidence for either case. 

Experimental support for the Batchelor-Goulian model 

To test the predictions of their model, Batchelor and Goulian 
measured changes in the transcriptional activity of OmpR- 
controUed genes using a two-fluorescent reporter strain, which 
provided indirect evidence for concentration robustness of OmpR- 
P. Recently, Gao and Stock direcdy confirmed the predictions of 
the Batchelor-Goulian model in the PhoR/PhoB system using a 
Phos-tag based method allowing for a quantification of the PhoB-P 
levels as a function of total PhoB amounts [26]. Experiments were 
performed with the wild-type (WT) system as well as with a PhoB 
mutant (PhoB^^"'') which exhibits reduced interaction strength 
(affinity) with PhoR. Both measurements could be well described 
by Eq. (1) with a CJCp ratio varying between 0.1-0.2 (Fig. 3A, 
solid lines). Overlaying the response curves with the respective Cp 
values (dotted lines) indicates that the PhoR/PhoB system operates 
in the regime Ci«Cp since the threshold concentration {Cp), 
beyond which PhoB-P becomes constant, is approximately equal 
to the value of that constant, as expected from Eq. (2). The 
observed shift of the threshold concentration in the mutant strain 
results from the reduced affinity of PhoB^^"'' which is associated 
with a larger value for Kp. Since Cp '^Kp, increasing Kp leads to an 
increased value of Cp so that the asymptotically constant 
phosphorylation level of PhoB — P « C™^ is reached at higher 
total PhoB concentrations, i.e. for total PhoB» lOpmol (Fig. 3A). 

Concentration robustness has also been observed in the 
reconstituted NRJI/NRI system of E. coli under in vitro conditions 
[27]. However, in that case the shape of the response curve is quite 
different (Fig. 3B): The dependence between [NRI-P] and total 
[NRI] does not appear to be linear below the threshold 



concentration and the asymptotically constant phosphorylation 
level («3.7|iM) is only reached for very large values of total 

[NRI] (»40nM). Together, this indicates that the NRII/NRI 
system operates in the regime Ci » Cp and, indeed, fitting the 
measurement data to Eq. (3) supports this view (Fig. 3B, solid line). 
Moreover, since in vivo concentrations of NRI are typically much 
lower than the threshold concentration of 40 |iM [28] it has been 
argued that, in the NRII/NRI system, concentration robustness 
will most likely not play a role under physiological conditions [27]. 

Results 

Ultrasensitivity in covalent modification cycles with a 
bifunctional enzyme 

To understand how ultrasensitivity may arise in TCSs with a 
bifunctional HK it will be helpful to analyze the consequences of 
reciprocal regulation in a related, but more simple system first. To 
this end, the reaction mechanism in Fig. 4A, which describes the 
reversible phosphorylation of a substrate 5 by a bifunctional 
enzyme E, is considered. The enzyme exhibits both, kinase (ii^) and 
phosphatase (Ep) activities, which catalyze the phosphorylation 
(S^S*) and dephosphorylatioii reactions {S*-^S), respectively. 
The transition between the two activity states is mediated through 
binding of an aUosteric effector L. For simplicity, it is assumed that 
Ej( has no phosphatase activity and, conversely, Ep has no kinase 
activity so that effector-binding effectively inhibits the enzyme's 
kinase activity and, concomitantly, activates its phosphatase activity. 
Note that this system is similar to TCSs with a bifunctional sensor 
kinase where the autophosphorylation and phosphotransfer reac- 
tions are replaced by a covalent modification (cf Figs. 2C and 4A). 
Also, the bifunctional converter enzyme is assumed to have just a 
single catalytic site, which is supposed to mimic the fact that the 
phosphotransferase and phosphatase activities of the sensor kinase 
in TCSs are also likely to occur on a single catalytic site [17]. 

The dynamics of this system is described by the set of ordinary 
differential equations (ODEs) 
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Figure 3. Experimental observations of concentration robustness in TCSs. Comparison between predictions of the Batchelor-Goulian 
model and measurements in the PhoR/PhoB [26] and NRII/NRI systems [27]. (A) Symbols denote measurements of PhoB-P as a function of total PhoB 
amounts in the wild-type system (open squares) and in a PhoB mutant strain (filled circles) (data were taken from Fig. 4C in Ref. [26]). Solid lines 
were calculated from Eq. (31) with CfT = I.2 pmol, C^~^ = 0.23 pmol and C™ = 10 pmol, C™" = 2.3 pmol. Note that Cp (dotted lines) determines 
both, the threshold amount of total PhoB beyond which PhoB-P becomes constant as well as the value of that constant, as expected from Eq. (2). (B) 
Symbols denote in vitro measurements of NRI-P as a function of total NRI (reproduced from Fig. 4A in Ref [27]). Solid line represents the best fit of the 
data to Eq. (3) with C, = 40 i^tM and Cp = 3J ftM, which indicates that the NRll/NRl system operates in the regime C,»Cp. 
doi:1 0.1 371/journal.pcbi.1 00361 4.g003 
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Figure 4. Ultrasensitivity in covalent modification cycles with a bifunctionai converter enzyme. (A) Reaction scheme: A substrate 
molecule (S) undergoes reversible phosphorylation by a bifunctionai converter enzyme which can exist In two activity states. Binding of the allosteric 
effector L Inhibits the kinase activity (Ek) and, concomitantly, activates the phosphatase activity (Ep) of the enzyme. (B) As the value of the 
dissociation constant K,i = kj/k^ Is lowered from Kj^l ftM»ET=0.1 /.iM to K,i = lO^^ ^iM«Et (from right to left) the steady state curve 
becomes ultrasensitive near the transition point L^ = 0.05 i^iM, as defined In Eq. (21). The solid lines were computed from the full model using Eqs. 
(4)-(7). Dashed lines were computed from the reduced models using Eq. (14) (right curve) and Eq. (18) (left curve). Parameters: Sj-^lO/iM, 
ki=k2=k^ =kT; = \/s, k+ =k+ =20/(^iM-s) so that if, = = 0. 1 ^iM, = l/(fiMi) and kj ^IQ-^s (for K,i = \0-^ ^iM) or /c,7 = l/.s (for 
K,t = lf,M). 

dol:1 0.1 371 /journal.pcbl.1 00361 4.g004 



^ = -^■+ [S*][Ep] [S*-Ep] +/c, [S-Ek] (4) 
=kt [S][Ek] - [K +h)[S-EK] 

^[^1=^ =ki [s*\m - {k^ +k2) [s- -Ep] 

^ = k+ [Ek] [L] - kj [Ep] - k+ [S*] [Ep] 

+ {k2+k2)[S*-Ep] 

together with the conservation relations for the total concentra- 
tions of substrate (Sp), converter enzyme (Ep) and allosteric 
effector [Lp) 



[S] + [S*] + [S- Ek] + [S* - Ep] = Sp 



[Ek] + [S- Ek] + [Ep] + [S* - Ep] = Ep 



(5) 



(6) 



d[Ep]/dt»0, which leads to the algebraic relations 



[S-Ek] 



[S][Ek] 



[S'-Ep] 



[S*][Ep 

K2 



[Ep 



[Ek][L] 



(9) 



Here, Ki=(ki+ A:f ) /k+ and K2 = (^2 + k^; ) /k+ denote Michae- 
lis-Menten constants associated with the kinase and phosphatase 
activities, respectively, whereas Kj = kj /k^ denotes the dissocia- 
tion constant for the enzyme-eflFector complex. 

Using the QSSA condition d[S* —Ep]/dtxO, it follows that 

d[S*] d[S*] d[S''-Ep] 



dt dt 
--h[S-EK]-k2[S*-Ep] 



(10) 



AEK][k 



-[S*] 



[S*][L] \ 
KiKd J 



where Eqs. (4), (8) and (9) have been used. In Eq. (10), [Ek] and [L] 
have to be found as functions of [S*] from the conservation 
relations (Eqs. 6 and 7) 



[L] + [Ep] + [S*-Ep]=Lp. 



(7) 



If the substrate concentration is much larger than that of the 
converter enzyme {Sp » Ep), one can neglect the concentrations of 
the enzyme-substrate complexes (since [S — Ek] + [S*—Ep]<Ep 
by Eq. 6) in the conservation relation for the substrate (Eq. 5), and 
the concentration of unmodified substrate can be expressed as 



(8) 



For later comparison, it will be useful to employ the quasi-steady 
state approximation (QSSA) in order to derive an effective equation 
for [S*]. By construction, the QSSA preserves the steady state 
structure of the underlying system [29] (which is the main focus 
here) although, for a better description of the transient dynamics, 
application of the total QSSA may be advantageous [30] . To apply 
the QSSA, it is assumed that, after a short transient period, the 
enzyme-substrate and the enzyme-effector complexes reach a quasi- 
steady state, defined by d[S — EK]/dtxO, d[S* —Ep]/dtxO and 



Et = [Ek]{1 



St-[S*] ^ [L] 
K\ Kj 



[S\ 



Et = [E][\ + 



[Ek] f, , [S* 



1 + 



K2 



(U) 



(12) 



Intuitively, it is clear that if the effector concentration is sufficiently 
large {Lp » Ep) the amount of effector that can be sequestered by 
the enzyme will be small since [Ep] + [S* —Ep]<Ep«Lp. Under 
this condition the conservation law for the effector (Eqs. 7 and 12) 
always reduces to Lp x [L] independent of whether the binding 
afifmity of the effector is high (if Kj is small) or low (if Kj is large). 
The latter only becomes important when the effector concentration 
is equal to or smaller than the enzyme concentration {Lp < Ep), e.g. 
under effector-limiting conditions. In the following, it will be shown 
that the type of effective equation, that is obtained from Eqs. (10)- 
(12), depends on the ratio Kj/Ep which may be regarded as a 
relative binding affinity for the enzyme-eflFector complex. 
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A low-affinity effector generates a graded response. If 

the relative binding affinity of the enzyme-effector complex is low 
{Et«Kii) one can neglect the terms associated with the enzyme- 
effector complexes in Eq. (12) since 



1 + 



< 



1 + 



«1 



provided that [S*]/K2 remains sufficiently small. Under this 
condition, one can use the simplified conservation relation 
Ltx[L] also at low effector concentrations (Lt<Et), so that 
[E]c] in Eq. (II) can be approximated by 



IEk\ 



1 + 



[SI 



Using this expression together with [L] xLj in Eq. (10) yields an 
effective equation for [S*] given by 



d[S*] 
dt 



St-[S'] 
Ki 



[S']Lt 
■ K2Kd 



5^ L./ [51 
Kd V K2 



(13) 



Under steady state conditions {d[S*\ldt = Q) the fraction of 
modified substrate exhibits a simple hyperbolic dependence on 
the effector concentration 



[S*] 
St 



Kd 



kiKi 
' kiK2 ■ 



(14) 



Hence, if reciprocal regulation occurs via a low affinity effector the 
system exhibits a graded response and ultrasensitivity cannot occur 
(Fig. 4B, right cun.'es). 

A high-affinity effector may lead to ultrasensitivity at low 
efifector concentrations. If the relative binding affinity of the 
enzyme-effector complex is high {Kd«Ej-) the simplified conser- 
vation relation Lt~[L] becomes invafid at low effector concen- 
trations {Lt<Et). In that case, the combination of Eqs. (11) and 
(12) leads to a quadratic equation for [Ek\, which can be written in 
the form 



1-^ 



1 + J 



i+j 



=0. (15) 



Here, x = [Ek] / Et and e = Kd / Et denote the normalized enzyme 
concentration and the relative binding affinity, respectively. In the 
limit e« I, one can neglect the 0(e) terms in Eq. (15) and obtains, 
to lowest order, the approximate solution 



[Ek] 



l + - 



(16) 



K, 



Note that this expression is only valid if the effector concentration 

is sufficiently small. The second branch of the solution (defined for 
Lt > Et) is of 0(e) and does not support ultrasensitivity (see SI 
Text SI). 



From the expression for [Ek] in Eq. (16) it follows that 
[Ex] /Kd~0{l/ e), i.e. [Ex] /Kd»l. Hence, one may approximate 
the free effector concentration (as determined by Eq. 12) through 



[L]: 



Kd 



K2 



(17) 



Finally, inserting Eqs. (16) and (17) into Eq. (10) yields the effective 
equation 



d[S*] 



ih(ET 



, St-[S*] 
Kx + St-[S* 



K2HS* 



(18) 



which is essentially the same equation as that derived by Goldbeter 
and Koshland for covalent modification cycles with two distinct 
converter enzymes [31]. Hence, if the activities of a bifunctional 
enzyme with a single catalytic site operate in saturation, so that 
rmx{K\,K2)<^ST, reciprocal regulation of the enzyme's activities 
by a high-affinity effector may result in zero-order ultrasensitivity 
similar as predicted by the Goldbeter-Koshland model (Fig. 4B, 
left curves). 

Comparison with the equation for the Goldbeter-Koshland 
model [31] 



rf[5* 
dt 



-kiKj 



St-[S*] 
Ki+St-[S*] 



-k2Pi 



K2 + [S* 



(19) 



shows that, in Eq. (18), the total kinase concentration (Kt) is 
replaced by Et — Lt whereas the total phosphatase concentration 
{Pt) is replaced by Lt. This result has an intuitive interpretation: 
If the binding affinity of the effector is sufficiendy high it can 
effectively sequester the enzyme into the states with phosphatase 
activity (\Ep] and [S*—Ep]) leaving only the enzyme fraction 
1 —Lt/Et for catalyzing the opposite reaction. In fact, using Eqs. 
(9), (16) and (17), it is straightforward to show that a high-affinity 
effector leads to a tight partition of the enzyme states according to 



p'PP : = [Ep] + [S* -Ep]kLt and 
K"f : = [Ek] + [S-Ek]kEt-Lt. 



(20) 



Hence, one may regard Lt and Et — Lt as apparent phosphatase 
and kinase concentrations, respectively. 

From Eqs. (16) and (18), it is also clear that ultrasensitivity 
becomes observable only at sufficiently low effector concentra- 
tions. Specifically, the transition from the 'on' ([5*]«5r) to the 
'off {[S*] «0) state, defined by k\{ET-L*j.) = k2L*j., happens at 



ki+k: 



(21) 



so that the occurrence of ultras(;nsiti\ ity is limited to the 
concentration range 0<Lt <Et, which is consistent with the 
range of validity of the approximation in Eq. (16). 

Batchelor-Goulian model with reciprocal regulation 

The Batchelor-Goulian model is based on the three activities of 
the sensor kinase shown in Fig. 2B, i.e. it essentially focuses on the 
signal transduction layer in the general scheme for two-component 
signaling depicted in Fig. 1 . However, within the context of this 
model it may become difficult to predict the input-output behavior 
as a function of the input stimulus, especially if the latter affects 
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multiple enzyme activities as observed in the PhoQ/PhoP/Mg^+ 
and NRII/NRI/PII systems (Fig. 2A). Guided by these examples 
the Batchelor-Goulian model will be extended by incorporating a 
mechanism that accounts for reciprocal regulation of the sensor 
kinase's autokinase and phosphatase activities by an allosteric 
effector. Analysis of this model shows that a low-aSinity effector 
may lead to stimulus-dependent concentration robustness whereas 
a high-afiinity effector may generate ultrasensitivity. In the latter 
case, the underlying mechanism is essentially the same as for 
covalent modification cycles (cf Fig. 4). 

To implement reciprocal regulation it is assumed (cf. Fig. 2C) 
that, in the absence of the effector, the free form of the sensor 
kinase {HK) can undergo autophosphorylation and mediates the 
phosphotransfer to the response regulator (step 1 and 2), but does 
not exhibit phosphatase activity (step 3). The latter is assumed to 
be activated through effector-binding (step 4), so that the 
phosphatase activity is carried by the ligand-bound form of the 
sensor kinase. Since HKl cannot undergo autophosphorylation 
(and phosphotransfer) binding of the ligand effectively leads to 
inhibition of the HK's autokinase activity and, concomitantly, 
activates its phosphatase activity. 

The dynamics of the extended model, as shown in Fig. 2C, is 
described by the five ODEs 



d\RR^ 
dt 



d[HKL] 
dt 



-- k, [RR - HK''] - k+ [RR''] [HKl] 
+k2[RR''-HKL] 



(22) 



that the ratio between the total concentrations of RR and HK is 
large {Rt/Ht» 1) [22,28] in which case one can use the simplified 
conservation relation (cf Eq. 8) 



[RR]xRt-[RR''] 



(30) 



instead of Eq. (27). Similar as in the case of covalent modification 
cycles (Eqs. 10-12), the steady state behavior of the system, 
described by Eqs. (22)-(30), depends on the affinity of the allosteric 
effector (Kj = kj /k^) relative to the total enzyme concentration 
{Ht). 

Note that for the derivation of Eqs. (22)-(30) it has been 
assumed that signal-sensing and the reactions describing the 
catalytic activities of the sensor kinase take place in the same 
compartment (the cytoplasm of the cell). Hence, this model 
directiy applies to cytosolic TCSs, such as the NRII/NRI system, 
but not to systems with a transmembrane sensor kinase where 
signal-sensing typically occurs in a different compartment. For 
example, in the PhoQ/PhoP system the sensor kinase PhoQ^ 
responds to changes of the Mg^+ concentration in the periplasm 
[20]. However, since effector-binding does not involve mass 
transfer the conditions for the occurrence of concentration 
robustness and ultrasensitivity are essentially the same (up to a 
factor accounting for the different compartment volumes) as those 
which are derived below on the basis of Eqs. (22)-(30) (see Methods). 

A low-a£Einity effector generates graded responses and 
stimulus-dependent concentration robustness. If the disso- 
ciation constant of the enzyme-effector complex is much larger 
than the total enzyme concentration [Kj»Ht) one can replace 
the conservation relation for the effector (Eq. 29) by [L] xLt, so 
that the steady state equation for [^^^ becomes (see Methods) 



(23) / K \ K 

ki[RR'^[HKL] + {k^+k,)[RR''-HKL] [RR'f -[RR^j^Rr + C^j^ + Qj +C,-^Rr = 0. (31) 



= ^k+ [HK] +k^ [HKP] +k, [RR-HK''] 



(24) 



-k^[L][HK]+k,[HKL] 



dt 



= k^ [RR''] [HKl] - {kj +kp) [RR'' -HKl] (26) 



together with the three conservation relations 

[RR] + [RR''] + [RR- HK''] + [RR'' - HK,] = Rt (27) 

[HK] + [HK''] + [HKl] + [RR-HK''] + [RR''-HKl] =Ht{28) 



However, this equation coincides with that, derived by Batchelor 
and Goulian in (Eq. 1), if the parameter Cp is substituted by the 
effective parameter 



(32) 



Hence, if the effector exhibits a low affinity ultrasensitivity cannot 
occur. Instead, Eq. (31) predicts a graded response of [^^^] with 
respect to changes in the effector concentration. 

To see this more explicitly, it wUl be useful to consider again the 
two hmiting cases Ci«C* and Ci»C*, which lead to the 
approximate solutions (cf Eqs. 2 and 3) 



C,«C;: [RR''] 



Rt 1 



C( \ Rt < Cp 



c;-Rt 



c;(i- 



c, 



Rt-c; Rt>C*p 



(33) 



[L] + [HKl] + [RR'' - HKl] =Lt 



(29) 



and 



where Rj, Hj and Lj denote the total concentrations of response 
regulator, histidine kinase and effector, respectively. Measure- 
ments in the PhoQ/PhoP and NRII/NRI systems have shown 



CIRt K*R 
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Figure 5. Stimulus-dependent concentration robustness in 
two-component systems. Steady state response curves according 
to Eq. (31) for C,«C; (A and B) and C,»C; (C and D). (A and C) [RR''] 
exhibits a graded response as a function of Lj. (B and D) [RR''] exhibits 
stimulus-dependent concentration robustness as a function of ^7-. The 
dotted lines indicate the threshold concentrations beyond which [RR''] 
becomes approximately constant. Note that, if ^7 = 1 mM (correspond- 
ing to the blue dotted line in A), increasing Rt beyond 4/iM does not 
lead to a higher phosphorylation level of the response regulator (B), 
which might explain why autoregulation in TCSs does not necessarily 
lead to a higher phosphorylation level of the response regulator (cf Ref. 
[22]). However, decreasing the effector concentration to LT = 03mM 
(corresponding to the red dotted line in A) allows [RR''] to increase as 
Rt increases. Solid lines were obtained from simulations of the full 
model (Eqs. 22-29) using the parameters: Rt^IO^M, Ht^OA jiM, 
k+ =k+ =2/UiM:s), k, = kp=k^ ^k^ = 1 /s, k^ =0.01 /(jiM-.s), 
kj =3/.s (.S:rf = 300^iM, cf Ref [46]). (A and B) k^ =0.1/*, /c+ =10/.? 
(C, = 0.1/i7W, Cp = 10/(M) and (C and D) kf^=50/s, k^^O.l/s 
(C( = 50 iiM, = 0.1 /.iM). Dashed lines correspond to the approximate 
solutions in Eq. (33) (A and B) and Eq. (34) (C and D). 
doi:1 0.1 371/journal.pcbi.1 00361 4.g005 



witli K*i = KdCpj(C,-\- Rt). From these expressions, it is apparent 
that [i^i^^] becomes constant at low effector concentrations and 
decreases as [i^i^^] ~ 1 /Lj if Lt becomes sufficiently large. IVlore 
precisely, if C, « C; (Eq. 33), [RR''] » Rt for Lt < C,,Kd /Rt and 
[RRP]'^CpKd/LT for Lt>CpKj/Rt (Fig. 5A). hi the opposite 
case, i.e. if Ci » C* (Eq. 34), the qualitative behavior of [i^i^^] is 
similar to that described by Eq. (33) although the transition from 
the state where [.Ri^^] is high (for Lt«K^) to the state where 
[i^ii^] is low (for Lt»K^) occurs more gradually (Fig. 5C). 

More importantly, concentration robustness is now predicted to 
occur in a stimulus-dependent manner since the maximal 
phosphorylation level of the RR {€*) depends on the effector 
concentration Lt (Eq. 32). However, since only C* (but not C,) is 
affected by Lt there is a notable difference between the two 
regimes, described by Eqs. (33) and (34), which may be used to 
distinguish them experimentally. In the first case, changing the 
effector concentration will change both, the threshold beyond 
which concentration robustness occurs and the value of the 
maximal phosphorylation level (both of which are determined by 
C*) (Fig. 5B). In contrast, when Ci»C* changing Lt only 



changes the maximal phosphorylation level while leaving the 
threshold concentration (which is determined by Q) unchanged 

(Fig.ro). 

Stimulus-dependent concentration robustness in the 
PhoQ,/PhoP system. Evidence for stimulus-dependent concen- 
tration robustness came from experiments with the PhoQ/PhoP 
system where IVIiyashiro and Goulian investigated the effect of 
genetic autoregulation on the expression level of PhoP-regulated 
genes at different ]V[g^+ concentrations in the growth medium 
[22]. At high ]V[g^+ concentrations, they observed almost no effect 
on PhoP-regulated genes indicating that the concentration of 
phosphorylated PhoP remained approximately constant under 
these conditions (despite an expected increase of the total PhoP 
concentration due to autoregulation of the /)/zoP gene). In contrast, 
under limiting ]VIg^+ concentrations a substantial increase in the 
transcript levels of PhoP-regulated genes was detected indicating 
that the PhoP-P concentration had increased under this condition. 

These findings were rationalized based on Eq. (1) by assuming 
that the parameter Cp, which determines the maximal phosphor- 
ylation level of the RR as well as the threshold concentration for 
reaching this level (Eq. 2), increases as the ]VIg^+ concentration 
decreases. Interestingly, such an inverse relationship between Cp 
and the effector concentration is readily predicted by the extended 
model (Eq. 32 and Fig. 5B), where it arises from the assumption 
that effector binding inhibits the autokinase activity and increases 
the phosphatase activity of the sensor kinase - in agreement with 
the observed regulatory effect of ]V[g-^+ on the activities of PhoQ. 
Indeed, in the opposite case, if effector binding activated the kinase 
and inhibited the phosphatase activity, Cp would be proportional 
to Lt- 

Although Miyashiro and Goulian did not measure the 
concentration of PhoP-P direcdy they observed a gradual (rather 
than switch-like) increase in the transcript levels of PhoP-regulated 
genes as the ]VIg^+ concentration was lowered - in qualitative 
agreement with the stimulus-response curves depicted in Figs. 5A 
and 5C. IVIoreover, measurements using isolated PhoQ sensor 
domains yielded an apparent dissociation constant for ]VIg^+ 
binding of x 300 iiM which is much larger than typical 
intracellular sensor kinase concentrations (i/7-<l|iM) [23,28]. 
Together, this supports the view that the PhoQ/PhoP system 
operates in the low-affinity regime {Kd»HT) described by Eqs. 
(33) and (34). Note that this conclusion is not affected by the 
circumstance that ]VIg^+ binding occurs in the periplasm. In that 
case, the low-affinity regime is characterized by the condition 
{K,i)^,»)'[Ht]^ (see Methods) where {Kd)^ denotes the dissociation 
constant of the enzyme-effector complex as measured in the 
periplasm, [Ht] ^ is the cytosolic concentration of the sensor kinase 
and y = Vcj Ve denotes the ratio between the cytosohc and the 
periplasmic volume. 

A high-affinity effector may generate ultrasensitivity at 
low efifector concentrations. If the dissociation constant of the 
enzyme-effector complex is much smaller than the total enzyme 
concentration [Kd « Ht) the steady state concentration of RR'' is 
determined by (see Methods) 



kfPP{HT 



Rt - [RR"" 



K^^^ + Rt-[RR''] 



ik„Li 



[RR'' 



KpHRR'']' 



(35) 



Here, the apparent catalytic rate of the phosphotransferase activity 
{k'l^'^) as well as the apparent Michaelis-Menten constant of that 
activity [K'^'^'^) are defined in terms of their intrinsic values [kt and 
Kt = {k,-\-k^)/k^) and the kinetic rates {k^ and k'^) associated 
with the autophosphorylation activity of the HK through 
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1 + 



k+ 



and KfP'^ = K,- 



(39) 



(36) 



Apparently, Eq. (35) is structurally identical to the steady state 
equation resulting from Eq. (18), so that the response of [-R-R^] 
with respect to Lj is predicted to become ultrasensitive if 



mz.x{Kt'^'^,Kp)«RT, 

and the transition from the 'on' state ([i^-R^] 
{[RRP]xQ) occurs at (cf Eq, 21) 



(37) 

; i^r) to the 'off state 
(38) 



Hence, for TCSs with a bifunctional sensor kinase the occurrence 
of ultrasensitivity is also restricted to low effector concentrations 
(0 < Lj < Ht) similar to covalent modification cycles with a 
bifunctional converter enzyme (cf. Eq. 21). However, compared to 
covalent modification cycles (Eqs. 18 and 19), the occurrence of 
ultrasensitivity in TCSs with a bifunctional HK only requires the 
phosphatase activity of the HK to operate in the zero-order regime 
(Kp«Rx). In contrast, the phosphotransferase activity can remain 
of first order as long as the regulatory factor, which multiplies Kf in 
Eq. (36), is sufficiently small, so that Kl'''' «Rt (Eq. 37). Thus, two 
scenarios are conceivable: First, if the regulatory factor is of order 
one or larger (/c^T > kt) both activities have to operate in saturation 
{rii2Lx(K,,Kp)<^RT) for ultrasensitivity to occur. Second, if the 
regulatory factor becomes sufficiently small, e.g. when 
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Figure 6. Ultrasensitivity does not require botli enzyme 
activities to be saturated. (A) As the phosphotransferase (PT) 
activity of the sensor l<inase changes from saturation (blue curve) to 
non-saturation (red curve) the steady state response of [RR''] as a 
function of Lt remains ultrasensitive, but the transition point (L^), as 
defined in Eq. (38), is shifted to lower effector concentrations. Blue 
curve: K,=0.1 fiM «Rj- = 10 j.tM, red curve: K, = Rt. In both cases 
isTf^^Cl/iAf. (B) Transient dynamics for Lr = 0.05 /iM (dotted line in 
A) indicating that the time-scale for reaching the steady state increases 
if the PT activity becomes non-saturated. Initial conditions: \RR\f^ = Rj, 
[HK^]g = HT, [L]g = LT, all other concentrations were set to zero. Solid 
lines were computed from the full model in Eqs. (22)-(29) with the 
parameters k+ =l.l/(fiM-s), ki = W/s, k^=OA/s, kj^=Om/s (red 



curve) and k+ =20/ (^iM-s), k,=ki 



--l/s (blue curve). Other 



parameters: Ht = \iiM, kj =20/(pM-s), kp=k^ =l/i, 
A:+ = IO/(/iM-.s), A:,7=0.I/.?, so that Kp = 0.\ tiM«RT and 
K,i = O.Ol iiM«Ht- Dashed lines in A correspond to the approximate 
expression for the stimulus-response curve in Eq. (40). 
doi:1 0.1 371/journal.pcbi.1 00361 4.g006 



Ki can become comparable to or larger than Rt without 
compromising the system's ability to generate ultrasensitivity 
(Fig. 6A). However, the condition on the kinetic rate constants in 
Eq. (39) leads to a shift in the transition point towards lower 
effector concentrations (Eq. 38) and may, substantially, affect the 
time scale on which the steady state is reached (Fig. 6B). 

Under the condition, stated in Eq. (37), the positive solution of 
Eq. (35) can be approximated by (see SI Text SI) 



RRn « 



Rt-K 



app 'V 



V 



K„ 



— Hf 



0<Lt<L*j- 



(40) 



which confirms that there is a sharp transition near the threshold 
concentration Lj, defined in Eq. 38, where the state of the 
response regxilar changes from almost full phosphorylation, i.e. 
[i^i^^Jwi^r for Lt<Lj, to a nearly unphosphorylated state 
[RRP]'^Kp«Rt for Lt>L*j. (cf Fig. 6A). 

Basal HK phosphatase activity may compromise 
ultrasensitivity. In the mechanism depicted in Fig. 2C it has 
been assumed that only the free form of the sensor kinase [HK] 
exhibits autokinase activity whereas ligand-binding has been 
required to activate the phosphatase activity (carried by HKi). 
However, it seems reasonable to also consider the more general 
case where the HK may exhibit some (low) phosphatase activity, 
even in the absence of effector. Conversely, HKl may also 
undergo autophosphorylation and mediate the phosphotransfer to 
some extent. To study the impact of such basal activities on the 
occurrence of ultrasensitivity it has been assumed that HK^ and 
HK catalyze the same set of reactions (Fig. 7A), but with lower or 
equal catalytic rate constants for the basal activities [kj<k,, 
kp < kp and k^ < k^). In general, changes in enzyme activity may 
also result from changes in the binding aflRnity. To account for 
such changes the association rate constants were allowed to vary 
according to k^ <k^ and k2 <k2 for enzyme-substrate 
binding and k^ f'-"' enzyme-effector binding (Fig. 7B). 

As can be seen in Fig. 7C increasing the basal autokinase and 
phosphotransferase activities of HKl, to the extent exhibited by 
HK, has only a minor effect on the response curve so that the 
occurrence of ultrasensitivity is not compromised in that case. In 
contrast, when increasing the basal phosphatase activity of HK 
ultrasensitivity gets lost if the affinity between HK and RR^ 
becomes sufficiently large (Fig. 7D). This suggests that, for 
ultrasensitivity to occur, RR^ must preferentially bind to HKf^ 
which requires tight regulation of the sensor kinase's phosphatase 
activity, e.g. through ligand-binding induced conformational 
changes of the sensor kinase [16]. 

Ultrasensitivity in the NRII/NRI system? Compared to 
PhoQ, which is a transmembrane sensor kinase, NRII is located in 
the cytosol where it controls the expression of nitrogen-regulated 
genes through reversible phosphorylation of NRI. The PII protein 
binds to the Idnase-domain of NRII which inhibits autophospho- 
rylation, but increases the phosphatase activity of NRII [16,21]. 
The components of the NRII/NRI/PII system have been purified 
and reconstituted in vitro [32] making this system amenable to 
measurements under well-defined conditions without interference 
from genetic autoregulation or other regulatory systems. 
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Figure 7. Impact of basal HK activities on the occurrence of ultrasensitivity. (A) Extended Batchelor-Goullan model (cf. Fig. 2C) with basal 
HK activities (dashed lines): Autophosphorylation (1), phosphotransfer (2 ) and phosphatase activity (3 ). (B) Scheme summarizing the allowed 
transitions between different enzyme states in the extended model with basal activities. It is assumed that ligand-binding occurs with reduced 
affinity (k^ <k^) if the sensor kinase has already undergone autophosphorylation (HK''). To prevent cycle fluxes under steady state conditions it is 
required that k'^ jk'^ I^k ■ Basal autokinase and phosphotransferase activities hardly affect the response curve. Even if the basal activities 
are identical to their full activities = ^j*" =/c,''"' =20/(|xM-.s), /cj" = 10/(nA/-.s), k^ =k'^' = l/s) the transition point (L^) remains 

almost the same and the steepness of the response curve is only slightly reduced (dashed curve). Basal phosphatase activity is assumed to be zero 
(k'^j = 0, k^ ' = 0)- The blue curve is the same as that shown in Fig. 6A where both basal activities are zero. (D) In contrast, upon increasing basal 
phosphatase activity the steepness of the response curve (ultrasensitivity) becomes substantially reduced. Dashed lines correspond to k'^ = 0.1kp and 
increasing values of A'j*^ =ak:^ for a = 0 (blue curve), 10^^, 10~^, 10~', 1. Basal autokinase and phosphotransferase activities are assumed to be zero 
(fc) =Q = k^' =ky ^^k')- Other parameter values are the same as for the blue curve in Fig. 6A. Simulations were done using Eqs. (64). 
doi:1 0.1 371/journal.pcbi.1 00361 4.g007 



In such a setting, Jiang et al. mea.sured the sensitivity of the 
steady state response of phospliorylated NRI (NRI-P) with respect 
to PII at different levels of total NRI (3-300nM) [27]. Half- 
maximal response occurred at [PII] 50 « 0. 1 |iM indicating that the 
Ki) for binding of PII to NRII is (much) smaller than the total 
enzyme concentration used in the experiments (total 
[NRII] = 1 |iAf). Hence, the two conditions Rt»Ht and 



Kd«HT, which are required for the apphcability of Eq. (35), 
seem to be fulfilled in the NRII/NRI/PII system, at least under in 
vitro conditions. However, even under saturating substrate levels 
(total [NRI] = 300 iiM) the response curve of NRI-P exhibited 
only a weak sensitivity with respect to changes in the PII 
concentration with an effective Hill coefficient of «j/wl.4 [27] 
(Fig. 8A). 




[PII], Lt [fJ,M] Lt [/uM] 

Figure 8. Autophosphatase activity of NRI may compromise ultrasensitivity in the NRII/NRI/PII system. (A) Comparison of experimental 
data (filled boxes, data taken from Fig. 4A of Ref. [27]) with the steady state response curve calculated from the extended Batchelor-Goulian model in 
Eqs. (22)-(29) with an extra term ' — k„p[RR'']' added to Eq. (22), which accounts for autodephosphorylation of NRI-P. The blue dashed line represents 
the best fit obtained for K, = OMiiM, Kp^l.84\iM, k,=OM/s and kp = 031/s. The other parameters were kept fixed: Rr^iOOuM, //r = l H^, 
k+ =k^ =k^ =kj =kj = l/s, k^ = 'lO/{iiM-s) so that K,i = 0.l nM and /i'„p = 0.0023/.? corresponding to a half-life of 5 minutes [27]. (B) As the 
autodephosphorylation rate constant of NRI-P is lowered (bottom to top: 0.0023/.?, 0.0005/.?, O.OOOl/i, 0/.?) the response curve becomes more and 
more ultrasensitive (solid lines). Note that ultrasensitivity is restricted to the region Lt<Ht as predicted by Eq. (38). The dashed (blue) lines in (A) 
and (B) are identical. 
doi:1 0.1 371/journal.pcbi.1 00361 4.g008 
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To explain this weak sensitivity Jiang et al. argued that the 
activities of NRII might operate in unfavorable kinetic regimes for 
ultrasensitivity to occur. Specifically, while the kinase/phospho- 
transferase activities were found to be saturated under the 
conditions of the experiments the phosphatase activity did not 
appear to be saturable, which is consistent with the theoretical 
prediction that ultrasensitivity requires saturation of the phospha- 
tase activity (c£ Eq. 37). However, given that intracellular NRI 
concentrations presumably lie in the submicromolar range [28] it 
seems unlikely that the condition KpKKRj- is still violated at total 
NRI concentrations a.s large as 300 [iM. This suggests that either 
conventional ideas about enzyme saturation are not applicable to 
the phosphatase activity of NRII (as discussed in Jiang et al. [27]) 
or that ultrasensitivity is compromised by another mechanism. 
The latter conclusion is supported by the observation (cf Fig. 8A) 
that the maximal phosphorylation level of NRI (~ 16 ^M) is much 
lower than the total NRI concentration used in the experiments 
(300 nM) which indicates the presence of a substantial phospha- 
tase activity, even in the absence of effector ([PII] <0.01 ^M). 

As indicated by Fig. 7D such an unregulated activity could 
result from a basal NRII phosphatase activity or, alternatively, 
from an intrinsic autophosphatase activity of NRI. Since the basal 
NRII phosphatase activity was found to be quite low [16] the 
second possibility appears more likely. In fact, compared with that 
of other response regulators the autophosphatase activity of NRI is 
comparably high [1] which results in a NRI-P half-life of 
5 minutes [27]. To study the impact of NRI-P autodepho- 
sphorylation on the occurrence of ultrasensitivity I have added an 
extra term { — kap[RR^]) to Eq. (22) and fitted the resulting set of 
equations to the measurements obtained by Jiang et al. under 
saturating conditions (Fig. 8A). To this end, only the Michaelis- 
Menten constants and the catalytic rate constants for (de- 
)phosphorylation were allowed to vary as these parameters should 
exhibit the most influence on the steady state response according 
to Eq. (35). The parameters Rt, Hj, and kgp were fLxed at 
their experimental values, whereas the remaining parameters (A:^ , 
, A:f , k2 and kj) were arbitrarily fixed at 1 /s so that they are 
all large compared to the autodephosphorylation rate constant 
kapX0.002'i/s. The thus obtained values for the Michaelis- 
Menten constants {Kp w 1 .84 nM and isTf « 1 .23 ilM) are much 
lower than the total NRI concentration (300 [iM) which suggests 
that the NRII/NRI/PII system operates in a kinetic regime that 
would, in principle, allow for ultrasensitivity. Hence, by lowering 
the autophosphatase activity of NRI the fitted response curve 
should become more and more ultrasensitive which is, indeed, 
what is observable in Fig. 8B. Together, this supports the view that 
the intrinsic autophosphatase activity of NRI might play a 
prominent role for the observed weak sensitivity of the NRII/ 
NRI system under in vitro conditions. 

Discussion 

In many two-component systems, the phosphorylation level of 
the response regulator protein is modified by a bifunctional sensor 
kinase which, apart from exhibiting autokinase and phosphotrans- 
ferase activity, also catalyzes the dephosphorylation of the response 
regulator through a phosphatase activity. In the present study, I 
have argued that the spectrum of potential input-output behaviors 
of such bifunctional systems does not only comprise graded 
responses [8-10] and concentration robustness [11,12], but also 
ultrasensitivity as it is well-known from phosphorylation-dephos- 
phorylation cycles with distinct converter enzymes [31]. To this 
end, I have proposed and analyzed an extension of the Batchelor- 
Goulian model [1 1] which considers the biologically motivated 



case where the autokinase and phosphatase ac:tivities of the sensor 
kinase are reciprocally regulated by an allosteric effector (Fig. 2). 

The analysis of the extended model showed that there exist two 
operating regimes under steady state conditions depending on the 
effector affinity: If the affinity is low compared to the total 
concentration of the sensor kinase {Kj»Ht) the system produces 
a graded response to changes in the effector concentration (Eqs. 33 
and 34) and exhibits stimulus-dependent concentration robustness, 
which means that the maximal phosphorylation level of the 
response regulator does not only depend on kinetic model 
parameters (as in tlu' original Batchelor-Goulian model), but also 
on the effector concentration. Consistent with experiments in the 
PhoQ/PhoP system [22], the extended model predicts an increase 
in the maximal phosphorylation level as the effector concentration 
is lowered (Eq. 32). However, if the effector affinity is sufficiendy 
high (K,i«Ht) the steady state equation for the extended model 
(Eq. 35) becomes structurally identical to that for covalent 
modification cycles with distinct converter enzymes (Eq. 19) so 
that ultrasensitivity may arise from the zero-order effect [31]. 

Apart from enzyme saturation due to the zero-order effect, 
sequestration of a signaling mokx-ulc into an inactive complex 
represents an alternative mechanism for the generation of 
ultrasensitivity in signal transduction networks [33-35]. Often, 
sequestration involves a reaction of the form 

'd k- 
A + B^AB, Kd=^ (41) 

where, by definition, A is sequestered by B into the complex AB. 
In this sense, regulation of enzyme activities by an allosteric 
effector may also be regarded as a form of sequestration. In the 
case of reciprocal regulation shown in Fig. 2C, the enzyme-effector 
complex (HKl) is not catalytically inactive, but rather has a 
different activity compared to the free form of the enzyme (UK). 
Buchler and Louis have shown that the simple mechanism in Eq. 
(41) can give rise to ultrasensitivity in the concentrations of A and 
B if the stoichiometric binding parameter Bt / Kj (where 
Bt = [B] + [AB]) exceeds unity, and the degree of ultrasensitivity 
increases as Bt/Kii»\ [36]. In the present study, the stoichio- 
metric binding parameter (Ht/Kj) plays a different role for the 
generation of ultrasensitivity since the condition Ht/Kii»1 does 
not guarantee the occurrence of ultrasensitivity per se, but only the 
validity of the reduced model, described by the steady state 
equation in Eq. (35). To obtain ultrasensitivity within the reduced 
model, the (apparent) Michaelis-Menten constants for the phos- 
photransferase and phosphatase activities of the sensor kinase also 
have to be sufficiendy small (Eq. 37), which distinguishes the 
mechanism, proposed in the present study, from purely seques- 
tration-based mechanisms. 

Interestingly, the idea of reciprocal regulation, as a mechanism 
to generate ultrasensitivity, does not seem to be restricted to two- 
component systems as the same mechanism may also apply to 
covalent modification cycles with a bifunctional converter enzyme 
(Fig. 4A). In both cases, reciprocal regulation may lead to 
ultrasensitivity if the stoichiometric binding parameters (Et / Kj 
in the case of covalent modifii:ation cycles or Hj- /Kj in the case of 
two-component systems) are sufficiently large. In this case, almost 
all free effector molecules are bound to the respective enzyme 
which leads to a tight partition of enzyme states into those with 
phosphatase activity and those with kinase activity (cf Eqs. 20 and 
52). As a consequence, the system behaves as if phosphorylation 
and dephosphorylation were catalyzed by independent enzyme 
subpopulations, which rationalizes why the corresponding steady 
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state equations (Eqs. 18 and 35) are structurally identical to that 
for covalent modification cycles with distinct converter enzymes 
(Eq. 19). However, this mechanism only 'works' as long as the 
enzyme is not saturated by the effector, which restricts the 
occurrence of ultrasensitivity to effector concentrations that are 
smaller than that of the respective enzyme (Figs. 4B and 6A). 

To assess the potential relevance of reciprocal regulation for the 
occurrence of ultrasensitivity under physiological conditions one has 
to evaluate to what extent the requirements for its occurrence 
(substrate excess, a large stoichiometric binding parameter and 
saturation of the sensor kinase's phosphatase activity) are satisfied in 
a particular system in vivo. Based on measurements in the EnvZ/ 
OmpR, PhoQ/PhoP and PhoR/PhoB systems, it seems that the 
requirement of substrate excess does not represent a limitation for 
the occurrence of ultrasensitivity as response regulator proteins are 
typically much more abundant than their respective sensor proteins 
[22,23,26]. In contrast, estimation of the stoichiometric binding 
parameter appears more difficult due to the limited knowledge on 
the range of input signals for a particular sensor kinase and their 
affinities relative to the total enzyme concentration. In general, 
histidine kinases may sense different signals (such as ions, 
metabolites, small peptides or auxiliary proteins) with widely 
different affinities [3]. Hence, it is conceivabie that the same system 
produces a graded response with respect to a low-aflBnity effector 
and an ultrasensitive response with respect to another effector with a 
high affinity. For example, apart from mediating adaptation to 
Mg^+ -limiting conditions the PhoQ/PhoP system is also involved 
in the regulation of bacterial viruience. This transcriptioncd 
program is initiated by antimicrobial peptides that seem to bind 
to the same periplasmic site in the sensor domain of PhoQ^as Mg^+ , 
but with a 1 00-fold higher affinity [3 7] , which could potentially shift 
the stoichiometric binding parameter into a regime where sigmoidal 
responses become possibie. 

The occurrence of uftrasensitivity aiso requires saturation of the 
sensor kinase's phosphatase activity which means that the 
Michaelis-Menten constant, associated with that activity, has to be 
smaller than the total concentration of the response regulator. 
Measurements in the EnvZ/OmpR system have shown that the 
dissociation constant for the EnvZ-OmpR complex is 5-fold smaller 
than the total OmpR concentration which indicates that enzyme 
saturation is, in principle, possible under physiological conditions 
[23]. However, the occurrence of ultrasensitivity can also be 
compromised by a sufficiently strong, unregulated phosphatase 
activity which may arise from a basal phosphatase activity of the 
sensor kinase (Fig. 7D) or from an autophosphatase activity of the 
response regulator. The latter might explain why the NRII/NRI/ 
PII system exhibits only a weak sensitivity with respect to changes in 
the effector (PII) concentration (Fig. 8B). Alternatively, it has been 
speculated that the observed weak sensitivity results from a non- 
saturable phosphatase activity of NRII [27] which is consistent with 
the prediction that ultrasensitivity requires the phosphatase activity 
to operate in the zero-order regime (Eq. 37). On the other hand, it 
has been shown that single mutations in the dimerization domain of 
a sensor kinase can substantially affect its interaction strength with 
cognate and even non-cognate response regulator proteins [26,38], 
which suggests that binding affinities between sensor kinases and 
response regulator proteins are highly evolvable. Hence, it is 
conceivable that one may employ directed evolution or site-directed 
mutagenesis to 'adjust' these binding affinities in a favorable range 
for ultrasensitivity to occur. In this sense, the results presented here 
may atso guide the design of synthetic regulatory circuits which aim 
to implement ultrasensitive response behavior at the level of two- 
component systems [39]. 



Methods 

Steady state analysis of Eqs. (22)-(29) 

Under steady state conditions, the right-hand sides of Eqs. (22)— 
(26) are set to zero so that summation of Eqs. (22) and (26) readily 
yields 



k, [RR - HK''\ = kp [RR'' - HKl] 



(42) 



Similarly, summation of Eqs. (23) and (26) leads to the steady state 
relation 



\HKr 



mm 



(43) 



where Kj = k^ /k'^ denotes the dissociation constant for the 
enzyme-effector complex. From Eqs. (24)-(26) together with Eqs. 
(30) and (43) one obtains the expressions 



-^^^J Kp K^^,^"^^ 



^f{[RR''])[HK'' 



where /([J^i{^]) is defined by 



kt 



1-f 



k, Rt-[RR'' 
k^ K, 



(45) 



whereas K, = (ki +ki)/ki and Kp = (k2 +kp)/k2 denote Mi- 
chaelis-Menten constants associated with the phosphotransferase 
and phosphatase activities of the sensor kinase, respectively. 

Using the expressions from Eqs. (43) and (44) in Eq. (42) and in 
the conservation relations, Eqs. (28) and (29), yields the set of 
algebraic equations 



k, Rt-[RR''] _ kp [RR""] [L]f^^ kt Rt-[RR'' 



kx Kp Kd 



and 



[HK"] I 



IRR'' 



[L](A+f([RR^])^'^ 



kx 



[RR'' 



K, 



[RR'' 



(46) 



(47) 



(48) 



from which the steady state concentrations [^^^], [HK''] and [L] 
have to be found. 
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Similar as in the case of a covalent modification cycle with a 
bifunctional enzyme the type of steady state solution, that is 
obtained from Eqs. (46)-(48), depends on the affinity of the 
allosteric effector. If this affinity is low [Kii » Hj) the concentration 
of free effector is approximately equal to the total effector 
concentration ([L] xLt). Replacing [L] by Lt in Eq. (46) readily 
yields the quadratic equation in Eq. (31) with Ct = (k^ /kt)Kf and 
C* defined in Eq. (32). 

In contrast, if the affinity of the effector is sufficiently high 
{Kii«Ht) the combination of Eqs. (47) and (48) yields a quadratic 
equation similar to that in Eq. (15) 



/ 



\ 



( l-Li/Hr 



K, 



(49) 



nmn) (i + ^^^^ +/( m )) (i + ^ 



where x=[HK^]/Ht and e = Kdl Hj denote the rescaled enzyme 
concentration and the relative binding affinity, respectively. In the 
limit e«l, the solution of Eq. (49) can be approximated by [40] 



[HK''] =xHt^'- 



, Lt<Ht. (50) 



With this approximation the concentration of free effector 
becomes (c£ Eq. 48) 



\L] = 



i+/([i^i^-])Miri 



K„ 



(51) 



according to (cf Eq. 20) 



[HKl] + [RR'' - HKl] xLt and 
[HK] + [HK''] + [RR - HK''] kHt- Ljs 



(52) 



so that Lt and Hj — Lt may be regarded as apparent 
phosphatase and kinase concentrations, respectively. 

Two-compartment model for regulation by an 
extracellular effector 

For TCSs with a transmembrane sensor kinase autophosphory- 
lation, phosphotransfer and dephosphorylation occur in the cytosol 
whereas signal-sensing typically takes place in the periplasm (for 
gram-negative bacteria) or directly in the extracellular space (Fig. 1). 
Hence, a proper model would have to distinguish at least 3 
compartments: The cytosol (where the response regulator is 
located), the plasma membrane (to which the sensor kinase is 
confined) and the extracellular space (where the effector is located). 
For gram-negative bacteria one would also have to consider a 
periplasmic compartment as many sensor kinases seem to respond 
to signals in the periplasmic rather than directly in the extracellular 
space [3]. Together, this makes it difficult to propose a generic 
model for TCSs that are regulated by non-cytosolic effectors which 
will, therefore, not be attempted here. 

Instead, to evaluate the impact of compartmentalization on the 
conditions for the occurrence of ultrasensitivity and concentration 
robustness it seems reasonable to consider (as a first approxima- 
tion) a simplified model where the reactions describing the 
catalytic activities of the sensor kinase occur in the cytosol (similar 
as assumed in the original Batchelor-Goulian model) whereas 
binding of the effector to the regulatory site of the sensor kinase 
occurs either in the periplasm or in the extracellular space. 
Because effector-binding does not involve mass transfer between 
the extracellular space (or the periplasm) and the cytoplasm the 
equations for such a two-compartment model are essentially the 
same as those for a single compartment (Eqs. 22-30) if the mass- 
balance equations are written in terms of average molecule 
numbers (rather than concentrations). The corresponding ODE 
system then reads 



Using this expression for [L\ in Eq. (46) yields the equation 
Rt-[RR'^ 



=k. 



Kp Ht—Lt 



Rt-[RR''\ k^( k, Rt-[RR''\ 

^+ K, ^kjy^k^ K, 



jn^P=k,n^_uj.p-k+nj^j^p-nHKL+k2 nnnP_uj^^ (53) 



dt 



nHK^ = kt riHK -ni-kj hhkj^ - k^ n,^p 
hhkl + {ki +kp)n^p_HKr 



dt 



hhk =—kx nHK + k^ n^f^p + k,njfj^_jjj.p 
- k^ uhk-il + kj hhkt^ 



which can be rewritten in the form shown in Eq. (35) of the main 
text. 

Similar to the case of covalent modification cycles it is 
straightforward to show (using Eqs. 43-45, 50 and 51) that a 
high-affinity effector leads to a partition of enzyme states 



jnj^j^_ijKP=kt nRR-nuj.p - (k^ +k,)n^_jjj.p 



Jf"RRP-HKL =k2 nj^p riHKT^ - (^2 + kp)nj^p 
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where nx denotes the average amount of species X (measured in 
mol). Compared to Eqs. (22)-(26) the second-order rate constants 
, k2 and kj are now measured in units of 1/ (mol-s), i.e. they 
are independent of the volume of the compartment in which the 
corresponding reaction occurs. In contrast, first order rate 
constants {kf , k^, k^, kj , kt and kp) have the same unit 
(1/s) as before. Mass conservation is now expressed in terms of 
molecule number conservation for the total amount of response 
regulator (urj.), sensor kinase {riHj.) and effector («Ly) as 



tlL + riHKL + "rrP -HKl="Lj-- 

Since the structure of Eqs. (53) and (54) is identical to that of Eqs. 
(22)-(29) it is clear that the conditions for the occurrence of 
concentration robustness and ultrasensitivity are identical in both 
cases if concentration-based quantities are replaced by their 
respective molar counterparts. 

Specifically, ultrasensitivity is predicted to occur if the amount 
of response regulator is much larger than that of the sensor kinase 
(nj{j.»«jyj.) and if the afiSnity of the effector is suHiciently high. 
The latter condition is now expressed as 



[LtI 



riLj. 
'1^ 



(59) 



is measured with respect to the extracellular (or periplasmic) 
volume Ve. In the case of an extracellular effector, one may think 
of Ve as the effective volume that is accessible to each cell in a 
population. In general, the effective volume decreases as the 
number of cells increases, e.g. due to cell growth. However, for the 
present purpose will be taken as a constant parameter. In 
addition, it is assumed that the extracellular space is a well-mixed 
compartment so that effector-diffusion can be neglected. 

Using the definitions in Eqs. (58) and (59), Eqs. (56) and (57) can 
be written in the form 



{Kp)^ + [RR''l 



(60) 



and 



(61) 



(55) 



where 



where the dissociation constant Kd ■ 



measured in mol. 



Under these conditions, the steady state amount of phosphorylated 
response regulator is determined by the analog of Eq. (35) 



kf^^{nHj.-nLj.) 



Kr+nR^-n^P^''''"^Kp+n^R 



"rrP 



(56) 



where k^^^ = k,/(l+k,/k+) and Kf^^ = K,{l+k^ /k+) 

/ ( \ Arkt/k'^') are defined by the same expressions as in Eq. (36). 
Similar as Kj, the Michaelis-Menten constants Ki = (kt + k^)/ki^ 
and Kp = (kp + k2 ) / are measured in units of wo/. Conversely, 
if the effector has a low affmity (Kd»nH-p) the steady state amount 
of RR^ is determined by the analog of Eq. (31) 



n^RRP - nRRP \nR^ + Cp^ + CA+ Cp^riR^ = 0, 



(57) 



where the rescaled Michaelis-Menten constants Cp = (k^ /kp^Kp 
and C( = (k^ lkt)Kt are defined by the same expressions as in Eq. 
(31). 

To analyze the impact of the compartment sizes on the input- 
output behavior one has to rewrite Eqs. (56) and (57) in terms of 
concentration-based rjuantities. For this purpose, the concentra- 
tions of the response regulator and that of the sensor kinase 



^rrP 



(58) 



are measured with respect to the cytosolic volume Vc, whereas the 
effector concentration 



y= - 



(62) 



denotes the ratio between the cytosolic volume and that of the 
extracellular (or periplasmic) space. Also, in Eqs. (60) and (61) the 
dissociation constant and the Michaelis-Menten constants have 
been rescaled according to 



mc=T' (^.)c=^' w.=§^. {cp)=^. 



(63) 



which gives them the conventional unit M = molll. The rescaling 
is motivated by the fact that, in a concentration-based description 
of chemical reactions, second-order rate constants have to be 
proportional to the volume of the compartment in which the 
corresponding reaction occurs [41], i.e. 
{k^)^=Vck^ and {k2)^=Vck2 giving them units of 
\/{M-s) = l /(mol-s). 

Similar as Eq. (35), Eq. (60) predicts that ultrasensitivity may 
occur at low effector concentrations {{Lt) ^<y[HT\^ if the afiBnity 
of the effector is suflBciently high {{Kj)^«y[HT\^. The latter 
condition follows from Eq. (55) using that Kd = Ve{Kd)g (Eq. 63) 
and Hhj^ = V^[Ht\^. (Eq. 58). Hence, depending on the volume 
ratio y the occurrence of ultrasensitivity may be favored (if y > 1) or 
suppressed (if y < 1) compared to a system that is regulated by a 
cytosolic effector (for which y=\). For example, if regulation 
occurs via a periplasmic effector y may vary between 1.5 and 4 
corresponding to a periplasmic volume fraction of 20-40% of the 
total cell volume [42]. In contrast, if regulation occurs via an 
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extracellular effector the volume ratio may be substantially smaller 
than 1 (y« 1) (especially at low cell densities) which would make 
the condition {Kj)^«y[HT]c less likely to hold and, therefore, 
suppress the occurrence of ultrasensitivity. 

Interestingly, Eq. (61) does not explicitly depend on the volume 
ratio. Hence, if reciprocal regulation occurs via a low-afEnity 
extracellular effector {{Kj)e»y[HT]^) the stimulus-response curves 
predicted by Eq. (61) are identical with those depicted in Fig. 5 if 
one replaces Kj and Lx by their extracellular (or periplasmic) 
counterparts {Kj)^ and respectively. 

Extended Batchelor-Goulian model with basal HK 
activities 

The response curves in Fig. 7C and 7D have been generated 
using the following set of equations (the corresponding reaction 
mechanism is shown in Fig. 7A and 7B) 



= k,[RR- HK''] +k',[RR- HKl] 

-k} [iLR^] [HKl] -kf [RR''] [HK] 
+ k2 ( [RR'' - HKl] + [RR'' - HK] ) 

[HK][L]-kJ [HKl] -k+ [RR''] [HKl] 
+k',[RR-HK[] 

+ (fcj- +ki,) [RR''-HKl] -k+'[HKL] +k^ [HKf] 



(64) 



d[RR''-HKL] 
dt 



--k+ [RR''] [HKl] - {k2 +kp) [RR''-HKl] 



dt 



kf [RR''] [HK] - [k^ +kp) [RR''-HK] 



^^[^^^ =kr[RR] [HKl] - (V +^;) [RR-HKl] 



where [RR-\, [HK''] and [L] have to be replaced using the conservation 
relations 



[RR] + [RR''] + [RR-HK''] + [RR-HKl] + [RR''-HK] 
+ [RR''-HKl]=Rt 

[HK] + [HK''] + [HKl] + [HKI] + [RR-HK''] + [RR-HKI] 
+ [RR''-HK] + [RR''-HKl] =Ht 

[L] + [HKl] + [HKl] + [RR-HKI] + [RR''-HKl] =Lt. 



^^ = -k+[HK]+k^[HK'']+k,[RR-HK'']-k^[HK][L] Supporting Information 

+ k^ [HKl] -kf [RR''] [HK] + ( k^ -h/t' ) [RR^-HK] ^ext SI Derivations and additional analysis. This file 

^ ' contains derivations of Eqs. (l)-(3) and (40) as well as the 

asymptotic analysis of Eq. (15). 
(PDF) 

=fc+' [HKl] -k^ [HKl] -kf [m M 
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